clear
hold on
y3 = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3651, 3601, 3550, 855, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];
x = 50:50:2000;
plot(x, y3, '--')%1

hold on
y2 = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3651, 3601, 2444, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];
plot(x, y2, '-^')%0.6

hold on
y1 = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3651, 2911, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];
plot(x, y1, '-x')%0.3

hold on
y4 = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3651, 3601, 3551, 3492, 577, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];
plot(x, y4, '-o')%1.5

hold on
y5 = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3651, 3601, 3551, 3501, 2898, 63, 0, 0, 0, 0, 0, 0, 0, 0, ];
plot(x, y5, '-d')%2

hold on
y6 = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3651, 3601, 3551, 3501, 3451, 3401, 2635, 0, 0, 0, 0, 0, 0, 0, ];
plot(x, y6, '-s')%10

hold on
y0 = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3636, 3506, 3373, 3216, 2973, 2166, 318, 0, 0, 0, 0, 0, 0, 0, ];
plot(x, y0, '-.')

plot(x,y0, '-.', x, y1, '-x',x, y2, '-^',x, y3, '--', x, y4, '-o',x,y5,'-d',x,y6,'-s');
legend('\beta = 0','\beta = 0.3','\beta = 0.6','\beta = 1.0','\beta = 1.5','\beta = 2.0','\beta = 10','Location','northeast','Orientation','vertical','Box','on');

xlabel('Attack size, k');
ylabel('Final system size');
box on;






